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Abstract 



Two-loop /3-function and anomalous dimension are calculated for N = 1 
supersymmetric quantum electrodynamics, regularized by higher derivatives in 
CN \ the minimal subtraction scheme. The result for two-loop contribution to the 

/^-function appears to be equal to 0, does not depend on the form of regularizing 
term and does not lead to anomaly puzzle. Two-loop anomalous dimension can 
be also made independent on parameters of higher derivative regularization by 
(D ' a special choice of subtraction scheme. 

43 ' 



1 Introduction. 



Investigation of quantum corrections in supersymmetric theories is a very important 
and complicated problem. In principle, supersymmetric theories have better ultravi- 
olet behavior, than nonsupersymmetric models. For example, in N = 2 Yang-Mills 
theory perturbative divergences are present only in one-loop diagrams. In principle 
it follows from the fact, that in supersymmetric theories axial anomaly and anomaly 
of energy-momentum tensor trace belong to the same supermultiplet |l|, 0, [| || . The 
axial anomaly is known to be completely denned by the one-loop approximation (5|. [| . 
while the trace anomaly is proportional to /5-function [[/J. Therefore, due to the super- 
symmetry the /9-function should be also defined by the one-loop approximation. The 
same arguments can be applied to AT = 1 supersymmetric theories. However, explicit 
calculations of radiative corrections show, that the /3-function in N = 1 supersym- 
metric models has contributions from higher loops || |9], |TD|, |TT|. This contradiction is 
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usually called "anomaly puzzle" and was investigated in a large number of papers, for 

example 0, H, HI II II 0, HI, HH HE H M M- 

Usually different proposals to solve anomaly puzzle require to fix the form of the 
/5-function in all orders of perturbation theory. For example, in theories with matter 
the /^-function should be related with the anomalous dimension. For the first time such 
/5-function was obtained by Novikov, Shifman, Vainshtein and Zakharov (NSVZ) from 
investigation of structure of instanton corrections f24jl . Later this result was checked by 



explicit calculations, which were usually made by the dimensional reduction technique 
f2~5fl . Two- loop /^-function, obtained in this regularization, is shown to coincide with 
a prediction, following from NSVZ exact expression. However, three-loop /3-function 
|26| , |27| , |28| does not agree with it. Nevertheless, the deviations can be removed by 
a redefinition of the coupling constant |29[, the possibility of such redefinition being 



highly nontrivial [|30| . In principle it is possible to relate DR scheme and NSVZ scheme 
order by order [|3T] in the perturbation theory. 

It is necessary to especially mention paper f2~3f , in which /5-function is shown to 
depend on the normalization of matter and gauge superfields. In particular, NSVZ 
/5-function can be obtained after a special rescaling of these superfields, which re- 
duces kinetic terms to the canonical normalization. Otherwise, (without rescalings) 
/3-function is argued to be completely defined by the one-loop approximation. So, it is 
really quite possible to obtain zero contributions of higher loops. The problem is how 
to calculate the corrections. In principle, it is possible to look for the regularization 
or renormalization scheme, in which /3-function is equal to the one-loop result or coin- 
cides with NSVZ exact /5-function. For example, the calculation of super Yang-Mills 



two- loop /3-function in the differential renormalization [32] was made in [33]. Another 
interesting possibility is using of higher covariant derivative regularization B5[. For 
the supersymmetric Yang-Mills theory the Lagrangian of the regularized theory was 



constructed in |36[ . For electrodynamics construction of the regularized Lagrangian is 
simpler, because instead of covariant derivatives it is necessary to use usual derivatives. 
However, calculation of diagrams, regularized by higher covariant derivatives is rather 
complicated. In particular, explicit calculation of the one-loop quantum correction [] 
for the (nonsupersymmetric) Yang-Mills theory was made rather recently [|37], [3£| and 
gives the same result as the dimensional regularization. In principle it is possible to 
prove, that one-loop calculations using higher covariant derivative regularization (cer- 
tainly, complemented by the additional regularization for one-loop diagrams) always 



give the same result as the dimensional regularization [3S]. Investigation of two- loop 
corrections in theories, regularized by higher derivatives has not yet been done. 

In this paper we try to understand features of higher derivative regularization in 
supersymmetric theories and calculate two-loop renormgroup functions for massless 
N = 1 supersymmetric quantum electrodynamics in this regularization using minimal 
subtraction scheme. 



^^Note, that introducing of a term with higher covariant derivatives does not lead to regularization 
of one-loop divergences. For the one-loop divergences it is necessary to use one more regularization, 
for example, introduce Pauli-Villars fields. 
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The paper is organized as follows: 

In Section H we introduce notations and remind some information about iV = 

1 supersymmetric electrodynamics. In the next Section ^ the considered model is 
regularized by adding of higher derivative term. After it we describe the quantization 
procedure for the constructed theory. Two-loop /9-function and anomalous dimension 
are calculated in Section |j. Agreement of the results with renormgroup equations is 
checked in Section || A brief summary and discussion are presented in Conclusion. 
Technical details of calculations, including expressions for all Feinman diagrams, can 
be found in the Appendix. 

2 Supersymmetric quantum electrodynamics. 

N — 1 supersymmetric massless electrodynamics in the superspace is described by 
the following action: 




Here <p an d 4> are chiral superfields, which in components can be written as 

<f>{y, 9) = <p(y) + 0(1 + l5 )iP(y) + U(l + l5 )9f(y); 

}{y, 9) = (p{y) + 9(1 + 75M2/) + \9(l + l5 )9f(y), (2) 

where y^ = + i9 7 ^ 7 ^9/2 are chiral coordinates, ip and (p are complex scalar fields, 
ip and if) are Maiorana spinors, which can be unified in a Dirac spinor 

* = ^((1+7b)^ + (1-75)^), (3) 

and / and / are auxiliary complex scalar fields. 

V is a real abelian superfield, which is a supersymmetric generalization of the gauge 
field. In components this superfield can be written as 

V(x,9) = C{x)+iV29^{x) + \{99)K{x) + ^9 l5 9)H(x) + h^l^A^x) + 

+V2(99)9(i lbX (x) + zfTsd^ixj) + - A (99)\D{x) - \d 2 C{x)). (4) 

The chiral superfield W a is a supersymmetric generalization of the field strength 
tensor and in the abelian case is defined as 

2 In our notations the metric tensor in the Minkowski space-time has the diagonal elements (1, -1, 
-1,-1). 
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(5) 




(6) 



Model ([!]) is invariant under supersymmetric gauge transformations 



V- V- l -{A + A+)- 





(7) 



where A is an arbitrary chiral superfield. In principle, it is possible to choose Wess- 
Zumino gauge, in which the superfield V is written as 



However, this gauge is not supersymmetric. That is why we do not use it for calculation 
of quantum corrections. 

3 Higher derivative regularization of N = 1 super- 
symmetric electrodynamics. 

To regularize model (jl]) by higher derivatives let us first modify its action by the 
following way: 



Note, that the considered model is abelian and the superfield W a is gauge invariant. 
Therefore, a regularizing term should contain usual derivatives instead of the covariant 
ones. 

It is convenient to introduce operators 





S 



S — Sq + S\ 





2 




1 



,2 



1 



D 



(10) 
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which satisfy identities 



D 2 D 2 + D 2 D 2 = -16II 1/2 <9 2 - 16<9 2 ; (11) 
~D^ 5 D DY15D + ^frflsD = -16rT n 1/2 d 2 - 160"0". (12) 

Then the first term in action (|9|) can be presented in the following form: 



1 „ aln 

^uge = ^Re J d*xd 2 6W a C ab (l + —)W b = 

1 r rP' n 

= -^J^^9VU 1/2 d 2 (l + —)v. (13) 



The gauge invariance (^) can be fixed by addition of the terms 

which are invariant under supersymmetry transformations. Then due to identity (|TT 
the kinetic term for the gauge field is written in the most simple form: 



(14) 



1 „ azn 

S gauge + S gi =—j d*x d*6Vd 2 (l + —)V. (15) 

Due to the gauge invariance (|7]) the renormalized action of the considered model 
(without gauge fixing term) can be presented as 



1 t & 

Sren = ^ 3 (A//i)Re J d A x d 2 6 W a C ah (l + — )W 6 + 

+iz(A//i) J d 4 x d 4 6 (fe 2V <\> + 4>*e- 2V 4>) , (16) 

where e = e(A//x) is a renormalized coupling constant. A bare coupling constant eo is 
defined by the equation 

3 = ^ZsCA/a*) (17) 

and does not depend on A///. The /^-function and anomalous dimension in our nota- 
tions are defined as 

_ d / e 2 \ rflnZ 
d In /i I 47r y ' (i In /x 

At the first sight, the generating functional can be written in the following form: 
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{1 c)^^ 
i — J d 4 xd A 6Vd 2 (l + —)V - 

-^(^(AM)-i) J d*xd*evu l/2 d\i + —)v + 



1 



+ ^Z(A//i) / d 4 xd 4 ^(0*e 2 V + 0*e-^0) + / d 4 xd 4 6JV + 



+ J d A x d 2 e (j (j) + j $) + J d A x d 2 e + 



3 



(19) 



(Below we slightly modify this expression.) Note, that the considered case corresponds 
to the gauge group U(l) and, therefore, diagrams with ghost loops are absent. 
Taking into account, that 



D 2 D 2 (f) = -169 2 
for any chiral superfield and that 

d 4 xd 2 9 = -- / de- 
generating functional ( |19"D can be presented as 



Z - / DV D<f> D<f> exp { i / d 4 x d 4 6 I -j^Vd 2 (l + — ) V 



(20) 



(21) 



1 pS,n i 

-— (z 3 (A//x) - l) VU 1/2 d 2 (l + ^)V+ -Z{k/^e 2V <P + ^e~ 2V $) + 

D 2 ~D 2 ~ D 2 ~ D 2 ~ \ 1 



8d 2J ^8d 2J T 8d 2 ' 



8d 2 ' 



(22) 

In order to calculate this functional we should present the argument of the exponent 
as a sum of a part Sq, quadratic in fields, and interaction Sj, where 

1 r ( f) 2n 1 

S Q = — j d*xd*9 [Vd 2 (l + ^)V+ ~(cf>*cf> + + 

D 2 ~ D 2 ~ D 2 ~ D 2 ~ \ 

00 1 t 

Si = E J d 4 xd 4 9 (<f>*(2VY<t> + ^{-2V) n $) - 

-(Z 8 (A/A*) -l)— 2 j d*xd*6VIl 1/2 d 2 (l + —)v + 

+ (Z{k/n) - l) i / ^ ^ (/ e 2V > + 0*e" 2y 0) . (23) 
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(Last two terms of Si will generate vertexes with insersions of counterterms, obtained 
in previous orders of perturbation theory.) 

Then the generating functional can be written as 



Z = exp^iSi(jjj, 7^---J| / DVD</>D4exp(iS Q ). (24) 

Note, that the differentiation over chiral superfields in our notations is defined as 
follows: 



6j(8 x ,x) D 



2 



5j(6 y ,y) 2 ~ ' 

so that 



-—ne x -9 y )8\x-y), (25) 
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d 4 xd z e x j(x,e x )<j ) (x,6 x )= ( j ) (y,6 y ). (26) 



The integral, remaining in equation (0), is Gaussian and can be easily calculated: 



z = exp {M\ij']Tr--)} x 



x exp [iJd^O (f± j + ~f±] - ^ 2(1+ ;1 /A2W ^ 



(27) 



Expansion of this expression in powers of J, j and j gives a series of perturbation 
theory. 

However, introducing of higher derivative term does not eliminate all divergences. 
Really, in Appendix |A| we check, that the superficial degree of divergence for model (|5[) 
is equal to 

uj a = 2 - 2n(L - 1) - E^n + 1), (28) 

where L is a number of loops and is a number of external 0-lines. Note, that u\ 
does not depend on a number of external Valines Ey. Therefore even after introducing 
of the higher derivative term with n > 2, divergences are present in one- loop diagrams. 
In order to regularize them || it is necessary to insert Pauli-Villars determinants in 
generating functional (|IDD, so that 



( ~ I 1 f d 

z = j dv n ( det pv ( v > M i)) 1 ex p r 4^2 y dAx die vd2 i 1 + J^) V 

__(z 3 (A/ A4 )-i) / d*xd±evu l/2 d\i + — )v + 
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+ 



^Z(A/fi) J d 4 xd 4 (0*e 2y + 0*e" 2y 0) + 

(29) 



• Jd 4 xd 4 9JV ■ j d'\rd' 2 (jo • jo) + / d'r d-0 +j*o 
where 



(det PV(V,M 



-l 



1 



£>$ L>$ exp Y 
1 



Z(A/n)^J d 4 xd 4 ($*e 2y $ + 



+ $*e~ 2y $ J + -Jd 4 x d 2 M$$ + - J d 4 x d 2 M$*$* 
and the coefficients q satisfy equations 

5> = i ; £qM 2 = o. 



(30) 



(31) 



Below we will also assume, that Mj = a* A, where a* are constants. Insersion of such 
Pauli-Villars determinants allows to cancel remaining divergences in all one-loop dia- 
grams, including diagrams with insersions of counterterms. 

Repeating the above arguments it is possible to write the Pauli-Villars determinants 
in the following form 



(detPy(V,M)) 1 = exp^i(S PV )i(v, ~,...^|exp|i j d 4 xd 4 



x 



x J 



1 



d 2 + M 2 



i 



tJ+J 



M D ~ 



d 2 + M 2 ° J d 2 + M 2 4d 2 



t .« M D 2 ^ 

'■' J a 2 + M 2 49 2j 



, (32) 



jj=o 



where 



00 1 r 

= J2—r d 4 xd 4 6 ($*(2F)"$ + $*(-2V) n <l) + 

+ (Z(A/fi) d 4 xd 4 ($*e 2y $ + $*e- 2V l»), 



(33) 



that allows to find their perturbative expansions. Then it is possible to calculate the 
generating functional Z according to the prescription 



'15 16 



1 • r O ' ' ' 



Z = n(detPV(iA M( )) C "exp{^ V!W , Uj 

( < / < ^ <rt ( > *^' +? ^- J 8>(l + 8»./A») 



x exp 



J 



(34) 
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The generating functional for connected Green functions in our notations is written 

as 



W = -ilnZ, (35) 
and the effective action is defined by making a Legender transformation 

r = w - j ' d^xd^ejv - J d 4 xd 2 6(]<p + ]4)) - J d 4 xd 2 e(]*<p* + j*0*), (36) 

where J, j and j should be expressed in terms of V, and through solving of the 
equations 

xr SW 5W ~ 8W 

V = JJ ] = 17 ; = 7T (37) 

Expressions for Feinman diagrams in the coordinate representation can be found 
expanding generating functional ([34|) and substituting the result into equations fl35]) 
and (|36|) . Certainly, after this procedure V will contain only lPI-diagrams. Expres- 
sions for Feinman diagrams in the momentum space can be then obtained by Fourier 
transformations. Performing the calculations we used this algorithm and tried to avoid 
direct application of Feinman rules in order to be complitely sure in the correctness of 
numerical factors for all diagrams. However, for the sake of completeness we formulate 
Feinman rules for the considered theory, which allow to verify structure of expressions 
for the diagrams. 

1. External lines correspond to a factor 

n / > / ■ (38) 

where the index E numerates external momentums. 

2. Each internal line of ^-superfield corresponds to 

2 '~ 5 4 (#i-# 2 ). (39) 



(fc 2 + zO)(l + (-l) n k 2n /A 2n 

3. Each internal line — 0* or — 0* corresponds to 

(Note, that in the considered theory the action is quadratic in matter superfields, that 
allows to formulate Feinman rules in a bit different manner, than for, say, Wess-Zumino 
model.) 

4. Pauli-Villars fields are present only in the closed loops. Each internal line $ — $* 
or $ — $* corresponds to 
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D 2 D I 5\e l -e 2 ). (4i) 



4(fc 2 - Mf + zO) 
Internal lines $ — $ and $* — $* corresponds to 

respectively. Also it is necessary to add — E Ci ^ or eacn c l° se d loop of Pauli-Villars 

i 

fields. 

t, i i , r dik 

5. Each loop gives integration over a loop momentum J ^ 4 . 

6. Each vertex gives integration over the corresponding 9: J d A 9. 

7. There are numerical factors, which can be calculated expanding generating 
functional ([34]). 

4 Calculation of two-loop renormgroup functions. 



Let us calculate two-loop /^-function and anomalous dimension for a model, de- 
scribed by action (||). In the two- loop approximation (3- function can be found after 
calculation of diagrams with Ey = 2, = 0, presented at Figure |I|. Note, that each 
graph at this figure corresponds to a diagram with internal 0-line, a diagram with 
internal 0-line and a set of diagrams with internal lines of Pauli-Villars fields. As an 
example in Appendix |B| we present detailed calculation of one-loop diagrams. Ex- 
pressions obtained for the other diagrams are collected in Appendix |C[ Each of these 
diagrams has the following structure: 



Terms proportional to J d 4 9V(—p,9)V(p,9) are not gauge invariant and should dis- 
appear after summing of all Feinman diagrams, that is very convenient for checking 
correctness of the calculations. The other terms can be written as 



-ReJd 2 9—^W a (-p,9)C ab W b (p,9) £ f^A). (44) 

' diagrams 

Having performed the calculations, in the Minkowsky space we obtained, that the result 
for two-loop contribution to the effective action, corresponding to the two-point Green 
function of the gauge field, can be written as 
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Alf = Re J d 2 6 ^iy a (p)^^(-p)(A_ loop + / 2 _ loop + /p V + /Konishi), (45) 



(27 

where (for simplicity we omit +i0 in propagators) 



iff d A k 1 f d A k 1 

/l-loop = -2{J (2tt) 4 k 2 (k + p y-^ Ci J (2<kY {k 2 _ M ?) ((k + v) 2 - ' ^ 



is a total one-loop contribution, including contributions of diagrams with internal lines 
of Pauli-Villars fields; 



/2-loop 



2 



d A k d A q (k + p + q) 2 + q 2 — k 2 — p 2 



(2tt) 4 (2vr) 4 k 2 (l + (-l) n k 2n /A 2n ) (k + q) 2 {k + p + q) 2 q 2 {q + p) 2 

(47) 

is a sum of diagrams ( f?2j ) - ( |97D with internal lines of and fields; 



/pV = e 2 ^Q J 



d A k d A q 1 

x 



X 



(2vr) 4 (2vr) 4 p(l + (-1)™A; 2 «/A 2 ™ 
(k + p + q) 2 + q 2 — k 2 — p 2 



k + qf- M 2 ) ((k + p + g) 2 - M 2 ) (g 2 - M 2 ) ((g + p) 2 - M\ 

AM 2 



+ 



(* + q) 2 - M 2 ) (g 2 - M 2 ) ((g + p) 2 - M 2 ) J 
is a contribution of diagrams fl92| ) - fl97|) with internal lines of Pauli-Villars fields; 



(4* 



ze 2 A v r d A k M 2 
/Konishi - " ln " 2- c > / (2^1 (P_M 2 ) 2 ((fc + P ) 2 -M 2 ) (49) 

is a total contribution of diagrams (H) - ([101]) with insersions of one-loop counterterms. 
(Mi = cijA are masses of Pauli-Villars fields.) 

Similarly, anomalous dimension can be found after calculation of diagrams with 
E v — 0, E<f, — 2, presented at Figure Q. Calculation of the one-loop diagram is 
described in Appendix |B[ Results for the other diagrams are presented in Appendix 0. 
The total two-loop contribution in the Minkowsky space can be written in the following 
form: 
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Arf = Jd A 8-^- (ffa 6) <p(- P , 6) + fr(p, 6) #(-p, ej) x 



(2k) 

d A k e 2 



x < i 



(^) A 2k 2 (k + p) 2 (l + (-l) n k 2n /A 2n ' 
d A k d A q e 4 



(2tt) 4 (2tt) 4 &V(A; + p) 2 (q + p) 2 (l + (-l) n k 2n / A 2n ) (l + (-1) V n M : 
d 4 fc d A q e 4 

(2tt) 4 (2tt) 4 fc 2 g 2 (g + p) 2 (k + q + p) 2 (l + (-l) n k 2n /A 2n ) (l + (-1) n q 2n /A 



+ 



ze 4 , A r d A k 1 /• d A k d A q 



X 



4tt 2 ft J (2vr) 4 p(A; + p )2(i + (-1)«A; 2 -/A2n) 7 (2tt) 4 (2tt) 4 

e 4 (A; + g + 2p) 2 

P(A; + p) 2 g 2 (g + p) 2 (A; + q + p) 2 (l + (-1)"& 2 »/A 2n ) (l + (-l) n q 2n /A 2n 
d A q e A ( f d A k 1 



( 27r ) 4 q 2 (q + p) 2 (l + (-l) n g 2n /A 2n ) (2tt) 4 P(£; + g) 



r d A k 1 \ 

" V ^ ■/ (^O 1 (F - M 2 ) {(k + q) 2 -M 2 ))' (50) 

Results (|45|) and (0) are evidently invariant under supersymmetry transformations, 
because they can be written as integrals from products of superfields (Q) and (|5|) over the 
superspace. The gauge invariance is absent in equation (0) because the calculations 
were made only for diagrams with Ey — 0, — 2. Adding of terms, corresponding to 
diagrams with arbitrary Ey and = 2, will certainly restore the gauge invariance. 

Divergent parts of the integrals in equations (|46|) - fl49l) are calculated in Appendix 
[B]. Using results, obtained there, one can conclude, that the sum of contributions ([46]) 
and (fl~7l) gives NSVZ result for the /3-function. Expression (|48D is shown in Appendix 
to be a finite constant and does not contribute to the /^-function. However, sum of 
diagrams with insersions of one-loop counterterms fl49|) is not zero and exactly cancels 
contribution (fT7|). Actually, equation ( |^ ) produces Konishi anomaly [i"0|, calculated by 



using Pauli-Villars regularization according to a method, described in According 
to this method, an anomaly is equal to a contribution of diagrams with internal lines of 
Pauli-Villars fields in the limit Mj — > 0, while contributions of diagrams with internal 
lines of usual fields are equal to 0. 

Divergent parts of integrals in equation ( pU|) are also calculated in Appendix [B|. 

Using results of Appendix [D] it is easy to verify, that counterterms, needed to cancel 



two-loop divergences in the minimal subtraction scheme can be written as 

1 A r f) 2n 
AS = - 16^ M a Re J ^ + A^) Wb + 
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d 4 xd 4 9 (0 

x { _ 

1 7T fX IT 2 fl IT 2 fl V ^ A 2 



a, A a 2 , 2 A a 2 , A/' , M 3\1 . . 



that corresponds to 

47T 2 7T , A 



-ln- + 0(a 2 ); (52) 

a 3 ). (53) 



e 2 , a(A//i) ' /x 

/ . \ , a . A a 2 . 2 A a 2 , A/^ , Mj 3\ 
Z(A//i) = 1 + - ln- + — In 2 -In- ^ Ci ln-^ + - +0 

V/ 7T (X TV 2 fl It 2 fX \i A 2J 



Therefore the two-loop /3-function and anomalous dimension of N — 1 supersymmetric 
quantum electrodynamics, regularized by higher derivatives, are written as 



2 

13 = — + 0{a% 

7T 

7 (a) = + ^ ( E ^ In ^ + ^ + 0(a 3 ). (54) 

In particular, two-loop contribution to the /9-function appears to be 0, so that the beta 
function is completely defined by the one-loop approximation. 

The anomalous dimension 7(a) in the two-loop approximation does not depend on 
n or, by other words, on a form of regularizing term. However, it depends on the ratios 
of Pauli-Villars masses to the constant A. Nevertheless, the dependence on Mj/A can 
be removed by addition of finite counterterms, proportional to In Mj/A: 



Srcn = ^Z 3 (A/ri Re J d 4 x d 2 6 W a C°»(l + —)w b - ^ £ c, In -± x 

x Re J d 4 x d 2 9 W a C ab W b + Z(A/fx) J J d 4 x d 4 6 [fe 2V + 4>*e~ 2V (f>) , (55) 

that corresponds to another subtraction scheme, in which anomalous dimension is equal 
to 

7(«) = -^ + ^ + 0(« 3 ) (56) 

and does not depend on both n and Mj/A. In principle, in this scheme it is possible 
to consider, that A — > 00, Mj/A — > 00 instead of Mj = a,A. 
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5 Comparing of the results with predictions of the 
renormalization group method. 



The obtained results can be checked by the renormalization group method. It 
is well known [42], that in renormalizable theories terms proportional to In /i/A are 
completely defined by one-loop counterterms. Therefore, it is possible to calculate such 
terms by renormgroup equations and compare them with the result of calculation of 
Feinman graphs. 

Using the notation 

t = In t (57) 
A 

for the considered model renormgroup equations can be written as 

Z(t) =exp{ J dt 1 (a(t))}; * = / T^T- ( 58 ) 
Because in the one-loop approximation the /3-function is equal to 

(3(a) = a 2 A + 0(a 3 ), (59) 

in the lowest order 

a{t) = a (l + Piaot + 0(ajj)), (60) 
where ao = a(0). Expanding the anomalous dimension in powers of a 

7 (a) = cryi + a 2 ^ 2 + 0(a 3 ) = 71 (a + ftajj t) + a 2 , 72 + 0(a$) (61) 
and substituting it to the first equation of (EST) , we obtain, that 



Z(t) = 1 + 7l a 1 + nfaaft t 2 /2 + l2 a 2 1 + 7l 2 a 2 1 2 /2 + O(ojj). (62) 

Taking into account, that according to the results of one-loop calculations 71 = 
and Pi = 1/tt, the function Z should take the following form: 



Z(A//i) = 1 + ^ In - - l2 a 2 In - + O(ajj) 

TV jj /i 

a , A a 2 , A , A 



l + -ln- + ^-ln 2 72 a 2 In- + 0(a 3 ). (63) 



7T /X 7T 2 fj, jJL 



Comparing this expression with equation (|53D we see, that terms proportional to 
ln 2 /x/A coincide, that can be considered as a check of performed calculations. 
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6 Conclusion. 



In this paper we calculated two-loop /^-function and anomalous dimension for N = 1 
supersymmetric massless quantum electrodynamics, regularized by higher derivatives. 
In particular, two-loop contribution to the /^-function is found to be and not to de- 
pend on the form of higher derivative term. As we mentioned above, this result follows 
from the fact, that the axial anomaly and the anomaly of energy-momentum tensor 
trace in the considered model belong to one supermultiplet. However, to obtain it we 
have to perform calculations using higher covariant derivative regularization. In princi- 
ple, this regularization (complemented by the Pauli-Villars regularization for one-loop 
diagrams) allows to perform easy calculation of diagrams with insertion of countert- 
erms, which are proportional to Konishi anomaly [[HJ and have nonzero contribution. 



Possibly the results of the paper allow to assume, that contributions of all higher loops 
to the /3-function of AT = 1 supersymmetric electrodynamics, regularized by higher 
derivatives, are also equal to 0. However, to be completely sure in it, it is necessary to 
calculate scheme dependent three-loop /^-function. 

Two-loop anomalous dimension is found to be independent on the form of higher 
derivative term if one renormalizes the coupling constant in this term. However, 7(a) 
depends on the ratios of Pauli-Villars masses to the constant A. Nevertheless, it is 
possible to make anomalous dimension completely independent on parameters of higher 
derivative regularization (n and M/A) if one introduces some finite counterterms in the 
renormalized action, that actually corresponds to a different choice of renormalization 
scheme. 

Note, that the result for /3-function does not contradict to NSVZ result 



a 2 / 

/3(a) = -3(1 -7(a)), (64) 

because in considered model /3-function depends on the normalization of the matter 
superfields ||23|| . In particular, after a scale transformation making matter superfields 



canonically normalized, it is possible to obtain NVSZ /3-function. However, in the 
present paper this statement was not checked by explicit calculations and we hope to 
make it later. 
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A Superficial degree of divergence for the super- 
symmetric quantum electrodynamics, regular- 
ized by higher derivatives. 



In order to calculate the superficial degree of divergence for an arbitrary diagram in 
the massless supersymmetric quantum electrodynamics let us introduce the following 
notations: 

L is a number of loops, 

ly is a number of internal IMines, 

1^ is a number of internal 0-lines, 

Ey is a number of external VMines, 
is a number of external 0-lines. 

We will start with the result for the superficial degree of divergence of massless 
supersymmetric electrodynamics without higher derivatives [^3] 



u = 2 - E+. (65) 

Adding of the higher derivative term changes only the propagator of V^-superfield, 
which will be proportional to k~ 2 ~~ 2n instead of k~ 2 in the usual supersymmetric elec- 
trodynamics. Therefore, in the regularized theory the superficial degree of divergence 
is equal to 

u A = uj — Inly. (66) 

Taking into account, that 0-lines are continuous, it is easy to prove the following 
identity: 



L = ly + 1 - (67) 
Therefore expression (|66|) can be finally rewritten as 

w A = 2-2n(L-l)-^(n + l). (68) 

In principle this result can be also obtaind from the Feinman rules in the superspace, 
but the derivation is more complicated. 

B Calculation of one-loop /3-function and anoma- 
lous dimension. 
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Expanding generating functional (p^) and substituting it to effective action (^) it 
is possible to find, that a one-loop diagram with Ey = and = 2 in the coordinate 
representation is written as 




^ J d 8 xid 8 x 2 0*(xi, 61) 0(x 2 , i 



r8 



d 2 (l + d 2n /A 2n ) 12 <9 2 



12- 



(69) 



where 



J d 8 x = J d 4 xd 4 6; 5 8 12 = 5\x 1 -x 2 )5\e 1 -e 2 ). (70) 
After Fourier transformation this expression in the Minkowsky space can be written as 



d *P X*, 



(27T) 



x 



fc 2 (l + (-1)«P™/A 2 ™) 



(27T) 



5 4 (0! - # 2 



1 



D 2 D 2 5\6 l - 2 ) 



/ /• , , „ d 4 p d 4 k 
2J ~ ' (2vr) 4 (2vr) 4 



„2 



A; 2 (l + (-l)«A; 2w /A 2ri )(A; + p) 5 



(71) 



where we used the identity 



5 4 (^ - 6 2 ) D 2 D 2 5\6 1 - 6 2 ) = U A {6 X - 6 2 ). (72) 
The integral over d 4 k in equation ([H]) can be calculated after the Wick rotation: 



d A k 



d 4 k 



2 J (2tt) 4 k 2 (l + (-l) n k 2n /A 2n )(k+p) 2 2 7 (27r) 4 A; 2 (l + & 2 ™/A 2 ")(A;+p) 2 ' 

(73) 

Then it is possible to perform calculation in the Eucliedian space and analitically 
continue the result for imaginary p°. 

The diagrams contributing to the one-loop /3-function can be considered similarly. 

From equations (p4[), (|35|) and (p6|) it is possible to find, that in the coordinate 
representation 




-i I d 8 x l V{x 1 ) 2 '^-5 8 ll = -iJ d 8 x l V(x 1 ) 2 ^5\x l -x 1 ) 



Ad 2 



(27T) 



v(p,e)y(-p,0) 



rf 4 A; 1 
J^k 2 ' 



(74) 



17 



where we take into account identity (|T2|). The corresponding diagram with Pauli-Villars 
fields in the coordinate representation is written as 




iYl°i J d s x 1 V(x 1 



n 2 n 2 

A(d 2 + M 2 ) 11 
d 4 p 



(75) 



Y:ciJd*6-^v(p,e)v(- P ,e) J 



d 4 k 



(2vr) 4 k 2 - M % 



2 ' 



(A bar at diagrams with Pauli-Villarse fields denotes a part of a vertex corresponding 
to $* or $*.) 

For the next diagram equations (Bl), (EB) and (KBT) give the following expression: 




DtD( D(D( 

4d 2 6 12 4 g 2 d 12 



■I y d 8 Xl d 8 x 2 V(x 2 ) ^5l 2 DlDHy{ Xl ) 



D 2 D 2 
4d 2 612 



d 8 x,d 8 x 2 V(x 2 ) ^5 8 2 (DlDl + [Dl D 2 ]) {VM^ 

D\D\ 



Ad 2 
1 

4d 2 ' 



2n2 

ti 8 

4d 2 6u 



d s Xl d 8 x 2 V(x 2 ) -^5 8 2 i D 2 l D 2 1 V(x 1 )^^8 8 2 + 



-dJiD ll ^D l V{x 1 )^^8l 2 - W{x x ) d^^&S 



d 2 



o 2 



12 



(76) 



Here we took into accout, that 



8\9 1 -9 2 )D ai ...D ak 5\9 1 -e 



for k = 0, 1, 2, 3 and used identities 



(77) 



[D 2 , D 2 



(7? 



AiD^D D %ll5 D D 2 = Aid^D 2 . 

Taking into accont equation (72) in the momentum representation expression (|76|) 
can be written as 



"16 J d ' 6 (2tt)* (2tt)4 k 2 (k + v) 2 V ^ 6) X 

x [D 2 D 2 + ADY l5 D k^ + 16(p + fc)%) V(-p, 9). (79) 

Let us note, that the integral 
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d 4 k \ 

(2tt) 4 k 2 (k + p) 
is proportional to p^. Therefore, it can be presented as 



I (2irYhHk + pY (80) 



T = PjJ^ T Vp, f d 4 k (k + pf -k 2 -p 2 _ 
M p 2 v 2p 2 J (2vr) 4 k 2 (k + p) 2 

1 r d A k 1 p^ r d 4 k ( 1 1 \ 

~ ~2 Ptl J (2n) 4 k 2 (k + p) 2 + 2f J J2^j*[k?~ (k + p) 2 )' (81) 

The last term can be omitted, because with the corresponding contribution of the 
diagram with Pauli-Villars fields it is proportional to 



r d 4 k (\ v 1 \ r d 4 k ( 1 v 1 \_ 

J (2n) 4 \k 2 ^ Ci k 2 -M 2 ) J (27r) 4 \ (k + p) 2 (k + p) 2 -Mf) ~ 

(82) 

where we take into account, that both integrals are convergent and it is possible to 
make in the second integral a substitution k+p — > k. Therefore, the considered diagram 
takes the following form: 



= ~tJ §r^W v(p ' e) {DlDl + 1 2 lDl Dl] + 16{k+ pfk ^ v( - p ' e) = 

= 5 / Wj-'kH^W^- 9 » ~ < fc + ^ ~ ^ *>■ < 83 > 

This diagram has two corresponding diagrams with the loop of Pauli-Villars fields, 
presented below. The first diagram is calculated similar to the massless case: 

= \ !>/ d*x,d 4 x 2 V( Xl ) V(x 2 ) ^flz/ n x 

D\D\ _i r d 4 p d 4 k 1 
4(^ + Mf) ^ 2 ^ Cl 7 (2tt)« (2tt)' (p - Mf) ((k + p) 2 - M 2 ) 

xV(p, 0)( - d 2 U 1/2 + (k + p) 2 + k 2 )v(-p, 9). (84) 

Taking into account, that signs in <p*V<f) and <p*V<f) vertexes are different, the second 
diagram can be written in the following form: 
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4 / d* Xl d*x 2 V{x x ) V(x 2 ) J^Sf 2 x 



d 2 + M 2 



x5 4 (6 1 -9 2 



5 i2 = -^E^/ rf4 ^ 2 



d A P d 4 k M 2 

V(p, G{) V(-p, 9 2 ) x 



(27r) 4 (27r) 4 " ^'"^ " v k 2 -M 2 

d A p d 4 k 



(k + p) 2 -M; 



DlD 2 l 5\6 l -6 2 ) = -iY J c l jd A 6 



(2tt) 4 (2tt) / 



x 



xv(p,0)v(- P ,0) 



Mi 



{k 2 -M 2 )((k + p) 2 -M 2 )' 
Collecting the results we see, that the sum of the one-loop diagrams is equal to 



55) 



E C * 



d A k 



(2vr) 2 \k 2 (k + P y 



(k 2 -M 2 )((k+p) 2 -M 2 



<4 

i Re I d2d j^w a (p,B)c ab w b (-p,e) x 



(27T) 2 

d 4 k 



(2n) 2 \k 2 {k + p)' 



E 1 



- > c 



(fc 2 -M2)((fc+p) 2 -M2 



56) 



Note, that all noninvariant terms, proportional to V^ 2 , disappeared, that can be con- 
sidered as a check of the correctness of the calculations. 

After Wich rotation in the Eucliedian space the last integral over d A k can be written 

as 



1 r d 4 k 



E 



2 J {2-n) 2 \k 2 {k + p) 2 V \k 2 + M 2 )((k + p) 2 + M 2 



(87) 



This integral is calculated below in Appendix [B]. In order to find the contribution to 
the effective action the result should be continued for imaginary p and multiplied by 

d A p 



Re / d 2 6— ^ Wa (p,6)C ab W b (-p,6). 



(2*) 

Two-loop Feinman diagrams are calculated in the similar way. However, the calcu- 
lations are much more complicated and we do not describe them in details. 
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C Diagrams, giving nontrivial contribution to the 
two-loop /3-function and anomalous dimension. 



Below we present expressions for all Feinman graphs, giving nontrivial contributions 
to the two-loop /3-function. Each of these graphs corresponds to a set of diagrams, 
which consists of a diagram with internal 0-line, a diagram with internal 0-line and 
diagrams with internal lines of Pauli-Villars fields. In order to find contributions to 
the effective action it is necesary to add the factor 



/ 



(27T)*' 



Note, that for simplicity of notations we also omit +i0 in propagators. 



(89) 




d A k 



Vd 2 U 1/2 V 

1 



k 2 (k + p) 2 k 2 



-V 2 



1 



-V 2 



1 



(k+p) 2 -M 2 



(k 2 - M 2 ) ((Jfe + p) 2 - M 2 ) k 2 - M\ 



(90) 




d 4 k J\ v 1 \ 

(2tt) 4 U 2 ^° l k 2 -M 2 y 



d 4 k d 4 q 



x 



vd 2 u 1/2 v„ A{k t p + q)2 -f~ p2 - v 2 



(2tt) 4 (2tt) 2 k 2 (l + (-l) n k 2n /A 2n ) 

2 



x 



-E 1 



(/c + g) 2 (/c + p + q) 2 q 2 (q + p) 2 
rf 4 /c <i 4 g e 2 



{k + g) 2 g 2 



x 



(2vr) 4 (2tt) 2 k 2 (l + {-l) n k 2n /A 2n 

4(k + p + q) 2 -k 2 -p 2 -2M 2 



Vd 2 U 1/2 V x 



((* + <?) 2 - M 2 ) ((k + p + q) 2 - M 2 ) (q 2 - M 2 ) ((q + p) 2 - M 2 ) 

-V 2 



((* + q? - M 2 ) (q 2 - M 2 ) 



(91) 



(92) 
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/ 



d A k d 4 q o 

■ V 



(27r) 4 (27r) 4 k 2 (l + (-l) n k 2n /A 2n ) 
2 +E 



X 



(93) 



(k + g) 2 (g + p) 2 Y ' ((& + g) 2 - M 2 ) ((g + p)2 - M 2 ) . 




d A k d 4 q 



(2tt) 4 (2tt) 4 k*(l + (-l)"A; 2 VA 2n 
2 2 



x 



V 



+ 



(k + q) 2 q 2 (q + p) 2 \q 2 (k + q) 2 (k + q) 2 (q+p) 2 



^ r d 4 k d 4 q 
Y l l (27r) 4 (27r) 4 ^( 



X 



Vd 2 U 1/2 V 



1 + (-1)™A; 2 "/A 2 «) 
2(g 2 + M 2 ) 



x 



((fc + g)2 - M 2 ) (g 2 - M 2 ) ((g + p) 2 - M 2 ) 

2 2 (g 2 + M 2 )(g 2 + (g + p) 2 )-8M| 



2 q 2 



{(k + q) 2 - M 2 ) (g 2 - M 2 ) ((q + p) 2 - M\ 



(94) 




d 4 k d A q 



- Vd 2 Ii 1/2 V 



(2tt) 4 (2tt) 4 k 2 (l + (-l)™A; 2 VA 2n ) 
4 4 4 



x 



(A; + q) 2 q 2 (q + p)< 



+ 



^ /■ <i 4 /c ci 4 g 



(2tt) 4 (2vr) 4 p(i + (-1)«A; 2 ™/A 2 «) 



g 2 (A; + g) 2 (fc + g) 2 (g + p) 2 / 
1/2I 

1 



Vd 2 Tlu<?V x 



x- 



+V 4 



((fc + g)2 _ M 2 ) (g 2 - M 2 ) ((g + p) 2 - M 2 ) 
4 4 \ 



+ 



(g 2 - M 2 ) ((* + g) 2 - M 2 ) ((fc + g) 2 - M 2 ) ((g + p) 2 - M 2 ) J \ 



(95) 




o? 4 /c d A q 



(27r) 4 (27r) 4 A; 2 (l + (-l)«A; 2 ™/A 



2ra^ 



\/ 2 



q 2 (k + q)' 



+ 
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k 2 {k + p) 2 / 
k 2 + M 2 



(q 2 - M 2 ) ((k + q) 2 - Mf) 



; (96) 



d 4 k d A q 



(27r) 4 (2vr) 4 k 2 (l + {-l) n k 2n /A 

2(q 2 + M 2 ) 



<7 2 (A; + g) 2 

; (97) 



(g 2 -M 2 ) ((A; + g) 2 -M 2 )J 



ie 2 , A 
m — 



d 4 k 



An 2 fiJ (2tt) 4 
2 A r d A k 



Vd 2 U 1/2 V 



k 2 (k + p) 2 



4tt 2 /W (2tt) 4 
V 



Vd 2 U 1/2 V x 



2 (A; 2 + M 2 ) (k 2 + (k + pf) - 4k 2 M 21 



{k 2 - M 2 ) 2 ((k + p) 2 - M 2 ) ' {k 2 - M 2 ) 2 ((k+p) 2 - M 2 ) 



(98) 




ie 2 , A /■ d 4 A; 
In 



/ 



4tt 2 /ii (2tt) 4 
2 A r d A k 



vd 2 u l/2 v- 



l 



- V 



/i7 (27T) 4 
1 



k 2 (k+ P ) 2 

Vd 2 U 1/2 V x 
1 



(P-M 2 )((A; + p) 2 -M 2 ) \k 2 -M 2 (k+p) 2 -M 2 , 




W2 \n±[^-V 2 i 



An 2 ii J (2tt) 4 



ie 



In 



A r d A k 



4tt 2 /ii (2tt) 



k 2 



P -M 2 



(99) 



jfc 2 Y (k 2 -M 2 ) 2 )' 



(100) 



(101) 



Expressions for diagrams, giving nontrivial contributions to the two-loop anomalous 
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dimension are presented below. Again, in order to obtain corresponding contributions 
to the effective action it is necessary to add the factor 




% r d A k 



(2tt) 4 p(l + (_l)njfc2n/A2n^ (jfe + p )S 



(103) 




d 4 £; rf 4 g 



(2tt) 4 (2tt) 4 g2 (l + (_l)n ? 2n/ A 2n^ ( g + p )I 



-E 



X 



+ 



+2 



/ 



(i 4 fc <i 4 g 



+ 



/ 



(27r) 4 (2tt) 4 + (_l)n g 2n/ A 2n) ( g + p )2 \k 2 

d A k d A q e 4 / 1 

(27T) 4 (2tt) 4 ^ ~ (_l)n g 2n /A 2n) 2 + if 



' -E^t^tto I + 



k 2 -M? 



(104) 



(A; 2 -M 2 )((A; + g) 2 -M 2 )/' 



= -2 



rf 4 A; d 4 g 



(27r) 4 (2tt) 4 ^ + (-I)«g2n/ A 2n) ( g + p )I 



X 



x 'i-E c 1 



/c 2 



k 2 -M?r 



(105) 



/■ d 4 A; rf 4 g 
~J (2 7 r) 4 (27r) 4 g4 ( 1 



-E< 



+ (-I)«g2n/ A 2n^ 2 V& 2 (& + ?) 2 
1 \ 



(P-M 2 )((A; + g) 2 -M 2 )/' 



(106) 



24 




J (2tt) 4 (2vr 



(2tt) 4 (2tt) z 



x 



(107) 



x 



k 2 (l + (-1)™A; 2 ™/A 2 ")g 2 (l + (-1) "g^/A 2 ") (g + p) 2 {k + q + p) 2 ' 




(108) 



<i 4 A; <i 4 g 



(2tt) 4 (2tt) 4 p(i + (-1)«A; 2 "/A 2 «)g 2 (l + (-1) «g 2n /A 2n ) (g + p) 2 (A; + g + p) 2 ' 




(109) 



rf 4 A; rf 4 g 



(2tt) 4 (2tt) 4 + (-1)"A; 2 ™/A 2 ™)g 2 (l + (-1) ™g 2ri /A 2n ) (<? + p) 2 (A; + p) 2 ' 




o? 4 /c <i 4 g 



x 



x 



(2tt) 4 (2tt) 4 
e 4 (g + k + 2pf 



(110) 



k 2 (l + (-1)™A; 2 ™/A 2 ™)g 2 (l + {-l) n q 2n / A 2n ) (k + q + p) 2 {q + p) 2 {k + pf 




i A r d A k 
In — 



8tt 2 //7 (2vr) 4 A; 2 (l + (-l)™A; 2 VA 2 ")(A; + p) 2 ' 

(111) 




i , A r d^k 
In 



8tt 2 ^7 (2tt) 4 A;4(i + (_1)"A; 2 VA 2 «) 



+ (H2) 



25 



% A f d 4 k 
H . In 



8tt 2 // J (2tt) 4 £;2(i + (-l)«A; 2 VA 2n )(A; + p) s 




i A r d 4 k 
In — 



4tt 2 fiJ (2n)*k 2 (l + {-l) n k 2n /A 2n ){k + p) 2 ' 



(113) 




i A r d 4 k 
In — 



8tt 2 /i7 (2vr) 4 A; 4 (l + {-l) n k 2n /A 2n 



(114) 



D Calculation of integrals, regularized by higher 
derivatives. 



Let us calculate divergent parts of integrals, which were encountered in equations 
(45) and (|50D. First it is necessary to perform the Wick rotation. After it the integrals, 
which are present in equations ([45|) and ([50]) will be proportional to 



h 
h 
h 
h 
h 
h 
h 



d 4 k- 



, 2n 



k 2 (k + p) 2 (l + k 2n /A 
d 4 k d 4 q 



(k 2 + M 2 )((k + p) 2 + M 2 



(k + p + q) 2 + q 2 — k 2 — p 2 



d 4 k 



k 2 (l + k 2n /A 2n ) (k + q) 2 (k +p + q) 2 q 2 (q + p)< 
M 2 



(k 2 + M 2 ) 2 ( (k + p) 2 + M 2 



d 4 kd 4 q 



d 4 k d 4 q 



d 4 k d 4 q 



k 2 {k + p) 2 q 2 {q + p) 2 (l + k 2n /A 2n ) (l + g2n/ A 2n 

(k + q + 2p) 2 

k 2 {k + p) 2 q 2 {q + p) 2 {k + q + p) 2 (l + k 2n / A 2n ) (l + g 2 "/A 2 «) 

1 

k 2 q 2 (q + p)\k + q + p) 2 (l + k 2n /A 2n ) (l + q 2n /A 2n 



26 



/ A ~7~~ rrWM); 

7 g 2 (g + p) 2 (l + g 2n /A 2n ) 



-/ 



d 4 Jfe d 4 g 1 

x 



(27r) 4 (27r) 4 fc 2^ + k ^/K^) 

(k + p + q) 2 + q 2 — k 2 — p 2 



x 



k + q) 2 + M 2 ) ({k + p + q) 2 + M 2 ) (g 2 + M 2 ) ((g + p) 2 + M 2 ) ' 



d A k d 4 q 1 

' U ' " ' (27T) 4 (27T) 4 jfc2(l + jfc 2n/ A 2n) X 

M 2 , x 

X 2 • ( 115 ) 

((fc + g) 2 + M 2 ) (g 2 + M 2 ) ((g + p) 2 + M 2 ) 

In order to calculate integral I\ it is possible to use four-dimensional spherical 
coordinates 



ki = k sin 6*3 sin 9 2 sin 6*i ; 
ki = k sin # 3 sin # 2 cos 6\ ; 
k 3 = k sin # 3 cos # 2 ; 
&4 = k cos 03. 

(116) 

and direct fourth axis along p^, so that the integrand will depend only on # 3 and 



OO 7T 7T 2w OO 7T 

y d 4 A; = J k 3 dk J d6 3 sin 2 6 3 J d6 2 sm6 2 J d6 1 = 4tt J k 3 dk J d6 3 sin 2 6 3 = 



oo 1 

= x = cos0 3 ] = 47r / k * dk J dxVl-x 2 . (117) 

-1 

Taking into account, that k^p^ = kp cos 9 3 = kpx, the integral can be written as 



d 4 k 



k 2 (k + p) 2 {l + k 2n /A 2n ) 

OO 1 

= An J kdk J dx 



VT^ 2 



o 

OO 



(k 2 + 2kpx + p 2 ) (l + k 2n /A 2n ) 

2tt 



[kdkidx a/1 f r , (118) 

(k 2 + 2kpx + p 2 )(l + k 2n /A 2n j 

27 



o c 



where the contour C is presented at Figure |3]. The integrand here has singularities at 
branch points x — ±1, a pole x = oo and a pole 



x 



k 2 + p 2 
2kp 



Then it is easy to see, that 



(119) 



dx 



c 



k 2 + 2kpx + p 2 



2iri Res 



vT 



ar 



-27ri Res 



x A 



k 2 + + p 2 
, x = x ) = 2iri 



x = oo 



/c 2 + 2fcpx + p 2 
Therefore, the integral over angles is reduced to 



k 2 + p 2 \k 2 — p 2 \ 
i . , - - + z- 



Ak 2 p 2 



Ak 2 p 2 



(120) 



k 2 + 2kpx + p 2 



7T 



7T 



P > k 



and finally 



P r k 1 7 1 1 
h = 2tx 2 \ dk — -, r + 2tt 2 / dk - -. 

J p2 ^ + A ,2n/ A 2n^) ./ fc ( 1 + A; 2 ' 



tT + o(l) H In 



7T 2 A 2 ™+p 2 ™ 



n p 2n ^ p ' 2' 

Integral I2 can be calculated using standard methods First, using an identity 



(l + k 2n /A 
2vr 2 (ln^ + l)+o(l). 



(121) 



(122) 



ab 



dy 



1 



(ay + 6(1 -y) 



2 1 



(123) 



it can be written as 



J d 4 k 



(k 2 + 2A;py + yp' 



E 1 



1 



- > c 



P + 2A;py + yp 2 + M| 



(124) 



Each of these integrals diverges, but their difference is finite. Therefore, to simplify 
calculations it is convenient to use an auxiliary regularization, for example, the dimen- 



sional regularization. Then the integrals in equation (|124|) can be easily taken: 
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I 2 = lim 7i dy 



} Jm r( 2 - D/2) ( 

r(2) ^ 



yO--y)F 



D/2-2 



D/2-2 



IT 



9 ^ / Mi 
2?r 2 V Cj In — + , 

T V p \ 



AM 2 
1 H ^-arctanh^ 



AM 2 + p 2 



(125) 



where we take into account, that = 1. 

i 

To calculate divergent part of the integral I 3 note, that 1% = J 3 (p/A) and due to 
the logarithmical divergence 



I 3 = ai In 2 — + a 2 In — + V h — . 

P P £S V A / 



If ai = 0, then it is possible to find 

dh 



CI2 = lim 



p^o din A 



d 4 k d 4 q 



An k 



2n~2 



Q 2 + Ma* 



A^fl + Jfe^/A 3 ") (* + ?)V 



(126) 



(127) 



If this limit does not exist, then a\ 7^ 0. The integral in the right hand side of equation 
(127) can be taken, using four- dimensional spherical coordinates: 



<? 2 + Mm a 1 j [ j (<? + foe) \/I - x 2 



^ 5 77 \a~a = 4tt dq dx 



(fc + g)V V W (k 2 + 2kqx + q 2 ) 2 

1 00 / ^ 1 2 

=- 2 \hI d 4 w + = F / fa ^ = F (128) 

so that 

a 2 = Ann 2 / rf 4 /c ? = 4vr 4 , (129) 

7 A 2 "(l + /c 2 "/A 2 ™) 

Therefore, from equation (|126|) we conclude, that 

h = 4vr 4 ln- + 0(l). (130) 
P 
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In order to calculate integral I4 let us note, that 

if^I P~ sr = f(p/M). (131) 

J (2vr) 4 (p + M 2 ) 2 ((k + p) 2 + M 2 ) 



Therefore, instead of taking the limit M — > 00 it is possible to take the limit p — > 0, 
so that 



A (FfsT +o(1) = y +o(1) - (132) 

Divergent part of integral J 5 can be also easily calculated, because 

J 5 = / 1 2 = 47r 4 fln 2 - + ln-^ +0(1). (133) 
V P PJ 

To find a divergent part of I 6 let us consider 
d 



lim A — (j 5 - I 6 ) 



limAA/^ A r i _ (fe + g + 2pn x 
p^o rfA y V (fc + o + ») 2 / 



k 2 (k + p) 2 q 2 (q + p) 2 (l + k 2n /A 2n )(l + q 2n /A 2n ) ~ 

= lim/^A ~? + 9)P " V x 
P-oy (/c + g + p) 2 

4 2n/ A 2n 

x ^—L 2 = 0. (134) 

/c 2 (/c + p) V(g + p) 2 (l + P"/A 2 ™) (l + q 2n /A 2n ) 

(It is important to note, that all integrals here are convergent.) Therefore, 

J 6 = J 2 + 0(1) = 4vr 4 f In 2 - + In -) + 0(1). (135) 

V P PJ 

A divergent part of I7 can be calculated similarly: 

lim A-^-(/ 5 -27, 

p^o aA v 



lim A 4- / d*kd A q(—^- 
p^o rfA y \ (k + 



p) 2 (k + q + p) 2 J 

1 



x 



k 2 q 2 (q + p) 2 (l + k 2n /A 2n ) (l + q ^/ A 2n 
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" iA J (k + p) 2 (k + q + p) 2 

x 



d_ 

p"o "dA 



x 



k 2 q 2 (q + p) 2 (l + k 2n /A 2n ) (l + q 2n / A 2n ) 



= lim A— / d A kd 4 q 
p^o dA 1 



g 2 + 2(k + p)q- (q + pf 
{k + p) 2 {k + q + p) 2 

x 



X 



lim A-- / d A kd 4 q 
p^o dA 1 



k 2 q 2 (q + p) 2 (l + k 2n /A 2n ) (l + q 2n / A 2n ) 
2kq - p 2 



(k+p) 2 (k + q + p) 2 

x 



X 



2kq . d 
A 



~ I d ' kd4q W(k + q) 2 ^dA 
= 8nJ d 4 qd*k 



k 2 q 2 {q + p) 2 (l + k 2n /A 2n ) (l + q 2n /A 2n ) 
1 



(l + k 2n / A 2n )(l + q 2n / A 2n 

kq q 2n /A 2n 

k 4 q 4 (k + q) 2 (l + fc 2n / A 2n) (l + g 2n / A 2n) 



, (136) 



To calculate this integral we again use four- dimensional spherical coordinates and direct 
fourth axis along g M . Then similar to the case, considered above, the integral over angles 
is reduced to 



i 



x\f 1 — x 2 f X\J 1 — X 2 

47r / dx — — 2ix w dx == 

J k 2 + 2kqx + q 2 J k 2 + 2kqx + q 2 
1 c 

2 -j-, / xy/i — x 2 \ 9 t-. ( xy/l - x 2 \ 

= 4 " lReS [ k 2 + 2kqx + q 2 ' X = °°J tReS [ k 2 + 2kqx + q 2 ' * = X ° ) = 

A 2,(_J_ . ^(k 2 + q 2 ) 2 i\k 2 - q 2 \{k 2 + q 2 ) \ 

I Akq + 8Jfe 3 g 3 Y 1 > 



so that 



2tt 



x 



/dx 
k 2 + 2kqx + q 2 



7r 2 q 



7T 



2 k 



k>q; 



q > k. 



(138) 



Therefore, 
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d 7 
lim A—(h ~ 2/7) = -lQnn 4 / do — 
p-o dA V 5 v ^ ( 1 



q 2n /A 2n 



yi+q 2n /A 2n 

(00 g 

00 



— 4n7r 4 / da; 



x" 



+ / d£; — j- 



(l + x 11 ) 2 J 1 + y- n 
v 7 



— Ann 4 [ dx 



g 3 (l + A; 2 ™/ A 2 ™ 

l/as 



(1 + x n ) 2 J l+y n 
v ' 



— 4n7r / da; 



„n-l 



(1 + x n ) 5 



-4tt 4 



(139) 



and finally 



I 7 = -I 2 + 2tt 4 In - + 0(1) = 2tt 4 f In 2 - + 2 In - ) + 0(1). (140) 
2 p \ P P J 



Using equation (|125| ) integral J 8 can be written as 



I 8 = 2n 2 J2ci J d A q- 



q 2 (q + p) 2 (l + q 2n /A 2n 



( '1 M ^ 
x In h a 

V ? \ 



4M 2 
1 H — arctanh^ 



4M 2 + q 2 



, (141) 



where Mj = a«A, cij being constants. To calculate the divergent part of this integral 
let us consider first an integral 



A- 



f(A/q)=I f (A/p 



(142) 



q 2 (q + p) 2 ' 

where / is a function. Differentiating If over In A and setting then p = 0, we obtain, 
that 



A 



dA 



p=0 



|^ A ^( A/8 ) = _|^- / ( A/ ,) = 

= -27r 2 | dg-/(A/g) =2rr 2 (/(oo)-/(0)). (143) 
So, if the values /(oo) and /(0) are finite, then 



/ / = 27r 2 (/(oo)-/(0))ln^ + O(l). 



(144) 
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If the function / is taken to be 



fU/l) = E 



2tt 2 



(l + q 2n /A 2n f\ 
then from equation ( |144j ) we conclude, that 



4M 2 + q 2 ' 



(145) 



1 H ^-arctanhi 



g 2 (g + p) 2 (l + g 2 "/A 2n ) 2 \ 9 2 " V 4M 2 + q 



4vr 4 y Ci ln- + 0(l) = 4vr 4 In - + 0(1). 
i V V 



(146) 



However, it is impossible to substitute in equation (|144|) 

2tt 2 



/(a/?) 



(l + g^/A 2 ™) i 



2 E Ci ln y 



(147) 



because /(oo) does not exist. Nevertheless, the function / can chosen in following 
form: 



d 



2tt 2 



f[A/q =A— - 
V 7 rfA V (l + q 2n /A 



2 E c * ln — 



2n 

ln 2 nq 2n /A 2n _ , Mi 
"3 E °i ln T" + 
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2tt 2 



(l + g 2n /A 2 ™ 



9 



(l + g 2 "/A 2 ™ 



2 ' 



so that /(O) = and /(oo) = 2-7T . Then from equation (|143|) we obtain, that 



:i48i 



and, therefore, 



-2E^ln^ = 4vr 4 ln- + 0(l) 

g 2 (g + p) 2 (l + g 2 «/A 2 M < 9 P 



(149) 



2vr 2 | A 



g 2 (g + p) 2 (l + q 2n /A 2n 



, Mi „ 4 , ,A /, Ax 
2 E c * m ^T = 27r 4 ln 2 - + 0(ln-J. 
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P 



P- 



(150) 



Then it is necessary to calculate logarithmical divergences. For this purpose we subtract 
from integral (|1 50[) terms, proportional In 2 A/p and differentiate the result over InA: 
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r a d 
lim A— — 

p^o dA 



2n 2 J d A q 



1 x ^ , M 4 2 A 
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g 2 (g + p) 2 (l + g 2 "/A 2 ^ < 9 P 
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27T 4 + 47r 4 ^can^. 



(151) 



From equations (|141 ), (|146[ ), ( |15U| ) and ( J151[ ) we see, that the divergent part of J 8 is 
equal to 

Is = 2vr 4 ( In 2 £ + 2 In £ ( £ In f + 2) ) + 0(1) . (152) 
In order to prove, that integrals Ig and iio are finite at A — > oo, first note, that 



J 9 = J 9 (p/A); 



'10 



Ao(p/A). 



(153) 



Therefore, it is necessary to prove, that J 9 (p = 0) and iio(p = 0) are finite constants. 
Let us set p = and make a substitution 



(154) 



Taking into account, that M = aA, where a is a finite constant, the considered integrals 
can be written as 



2(K + Q)^Q^ 



'10 



d 4 K d 4 Q 

(27r) 4 (2VT) 4 #2^ + R 2n^ ^ K + Qy + a 2y(Q2 

d 4 K d 4 Q a 2 



2 i 



(27r) 4 (27r) 4 #2^ + K 2n^ ^ K + Q^2 + fl2 ) (Q2 + 



(155) 

(156) 



The integrals over Q are evidently convergent. However, it is necessary to check, that, 
after taking these integrals, the remaining integration over K will be also convergent. 
Possible divergences can arise at K — » or at K — > oo. In the limit K — > 
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d*Q 2{K + Q) li Q li f d 4 Q 2Q 2 1 



( 27r ) 4 ((i,+Q) 2 + a 2 ) 2 (Q 2 + a 2 ) 2 ■/ (2vr)MQ 2 + a 2 ) 4 24vr 2 a 2 ' 
rf 4 Q a 2 /■ rf 4 Q a 2 1 



((^ + Q)2 + a 2 )(Q 2 + a 2 ) 3 7 (Q 2 + a 2 ) 4 967r 2 a 2 ' 



(158) 



It means that, in equations ( |155| ) and (|156|) the integration over is convergent if 
X — ► 0. Similarly, in the limit K — ► oo 
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rf 4 g 
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{Q 2 + a 2 ) 3 32vr 2 (K 2 + a 2 ) 32vr 2 if 2 ' 



where we used equations ( |128| ) and ( |121| ). Therefore, due to the presence of higher 
derivative term the integration over K in equations ( 155| ) and (|156|) is also convergent 
at K — > oo. Thus the integrals Ig and I\q are proven to be finite in the limit A — > oo. 

Collecting the above results we can finally write integrals ( 115|) in the following 
form: 
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/ 8 = 2^(ln'£ + 21n£(£«lnf + ?)) +0(1); 
J 9 = 0(1); 

/io = 0(l). (161) 
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Figure 1: Feinman graphs, giving nontrivial contributions to the two- loop /5-function 
of N — 1 supersymmetric electrodynamics. 
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Figure 2: Feinman graphs, giving nontrivial contributions to the two-loop anomalous 
dimension of N = 1 supersymmetric electrodynamics. 
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Figure 3: Contour C for calculation of integral over x. 
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